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A pure S'-wave covariant model for the nucleon 
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Using the manifestly covariant spectator theory, and modehng the nucleon as a system of three 
constituent quarks with their own electromagnetic structure, we show that all four nucleon electro- 
magnetic form factors can be very well described by a manifestly covariant nucleon wave function 
with zero orbital angular momentum. Since the concept of wave function depends on the formalism, 
the conclusions of light-cone theory requiring nonzero angular momentum components are not in- 
consistent with our results. We also show that our model gives a qualitatively correct description of 
deep inelastic scattering, unifying the phenomenology at high and low momentum transfer. Finally 
we review two different definitions of nuclear shape and show that the nucleon is spherical in this 
model, regardless of how shape is defined. 



I. INTRODUCTION 

The elastic electron-proton polarization transfer exper- 
iments undertaken at the Thomas Jefferson National Ac- 
celerator Facility (JLab) P, d, ^ disclosed that the ratio 
of the electric Gsp to magnetic Gmp form factors of the 
proton is not constant as Q^, the square of the momen- 
tum transfer, varies (referred to as a lack of scaling). 
These were followed_up by the super Rosenbluth mea- 
surements at JLab [4] confirming older Stanford Linear 
Accelerator Center (SLAG) measurements (reanalyzed 
by Arrington Q) which originally showed that Gep and 
Gmp do scale. 

It now seems clear that the discrepancy between the 
JLab polarization transfer results and other measure- 
ments using Rosenbluth separation are essentially due 
to two-photon processes [1, 0, B Q ■ There are still some 
uncertainties about the effective sizes of two-photon ex- 
change effects but there is no doubt that the form factors 
extracted from the Rosenbluth method require significant 
corrections, while those extracted from the polarization 
transfer method require only small corrections. As a con- 
sequence the polarization transfer method is a more ac- 
curate way to determine the form factors, removing the 
possibility that the lack of scaling is an experimental ar- 
tifact. 

The lack of scaling was a surprise. While it had been 
predicted as long ago as 1972 [10] these predictions had 
been largely forgotten. The absence of scaling was seen 
by some as proof that the quark wave function of the pro- 
ton must have orbital angular momentum components 
L > 0, and this observation was further supported by 
arguments, based on light-cone wave functions, that the 
Pauli form factor, F2, must be zero unless L > compo- 
nents exist [ll|. In this paper we will discuss how these 
results depend on the light-cone formalism and show, us- 
ing the covariant spectator formalism, that (i) it is pos- 



sible to construct a pure S'-wave covariant wave function 
for the nucleon, and that (ii) L > components, while 
they may be a part of any realistic wave function, are not 
required to explain the data. We cannot conclude that 
either formalism is wrong, only that the the concept of 
the "wave function" is different. All of this is discussed 
in Sec. HVDl 

The JLab data has also stimulated discussion about 
the shape of the nucleon. Is the nucleon spherical or 
it is deformed? If distorted, is this a relativistic effect? 
The answers to these questions depend in part on how 
we define "shape" and will be discussed in some detail in 
Sec lylTBl 

Finally, in order to fully examine the implications of 
our covariant approach, we use it to calculate the quark 
distribution functions in deep inelastic scattering (DIS). 
We obtain reasonable, qualitatively correct, results. 

To set the stage for the physical approximations used 
in this calculation, consider Fig. [T] In the Nc ^ 00 limit 
(where Nc is the number of quark colors), crossed dia- 
grams are suppressed and the gluon exchanges between 
the qq pair interacting with the photon can be included 
either as part of a constituent quark electromagnetic form 
factor (usually described by vector dominance) or as con- 
tributions from higher Fock components of the nuclear 
wave function. (Since Nc — 3, corrections to this sim- 
ple picture are expected to be of the order of only about 
10%, and are probably small enough to be included when 
the parameters of our phenomenological model are ad- 
justed to fit the data.) At low energies and momenta, 
the description of Fig. [Ija) has the advantage that un- 
known short-range physics can be included in a few con- 
stituent quark parameters (such as the anomalous mo- 
ment). In this paper we will adopt the view represented 
by Fig. [IJa); our constituent quark form factors include 
the physics that comes from higher Fock components, ex- 
pected to be important at modest Q^, vanishing only in 
the DIS limit. 
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FIG. 1: (Color on line) Two views of the nucleon form factor. Left panel (a): the form factor is built from a constituent quark form factor 
and nuclear wave functions. Right panel (b): the quark is point-like, with its "structure" described by higher Fock states of the nuclear 
wave function; in this case the Fock state with 3-quarks and one gluon. 



The remainder of this paper is divided into five sec- 
tions. In the next section we define the nucleon wave 
function in some detail. In Sec. Ill, we derive the nucleon 
form factors using this wave function. The principal re- 
sults are Eqs. ([28| and ([29| . Readers uninterested in the 
details may skip directly to Sec. IV, where the first qual- 
itative results that can be derived from these equations 
are discussed. Then, before studying the form factors 
quantitatively, we use this model to derive, in Sec. V, 
formulas for the quark distribution amplitudes measured 
in deep inelastic scattering. In Sec. VI we parameterize 
both the nucleon wave function and the quark form fac- 
tors and determine the parameters by fitting the nucleon 
form factors and the moments of the quark distribution 
amplitudes. Finally, in Sec. VII, we review the work and 
draw conclusions. The appendices provide a few details 
considered too technical to be included in the text. 



II. RELATIVISTIC NUCLEON WAVE 
FUNCTION 

The relativistic nucleon wave function presented here is 
similar to the one developed by Gross and Agbakpe , 
referred to here as Ref. I. Two observations motivated 
us to develop this new 5- wave model. First, the overall 
motion of the composite nucleon described in Ref. I does 
not satisfy the Dirac equation for a spin-1/2 particle (we 
are not sure that this is an essential requirement, but it 
is certainly an advantage). Attempts to find a wave func- 
tion that does satisfy the Dirac equation lead us directly 
to this model. Second, as emphasized by Kvinikhidze 
and Miller [l^ , the model of Ref. I depends on the direc- 
tion of the relative three-momentum k between a pair of 
noninteracting quarks (referred to as a diquark although 
this term is usually applied to an interacting pair) and 
the third quark, and hence includes angular momentum 
components. The S'-wave model presented here is inde- 
pendent of the direction of this momentum, and exactly 
reproduces the pure S'-wave structure of the simplest non- 



relativistic SU{2) x SU{2) wave function when the nu- 
cleon is at rest. Since one of the goals of this paper is to 
study the implications of the absence of L > compo- 
nents in the wave function, it is essential that our model 
wave function have no such components. 

Following Ref. I , we consider the nucleon described 
by a wave function built in the spirit of the covariant 
spectator formalism [ij, [l^, which has already been 
solved and shown to work successfully for the three- 
nucleon system [l^ [l^, The nucleon with four- 

momentum P and mass M is described by a wave func- 
tion for an off-shell quark and an on-mass-shell diquark 
like cluster. 



*(P,/c) 



^i)-i(fc|r|p) 



(1) 



where F is the vertex function describing the coupling of 
an incoming on-shell nucleon with mass AI to an outgo- 
ing off-shell quark and an on-shell quark pair (the "di- 
quark"). The continuous mass of the diquark pair is fixed 
at some mean value, which scales out of the final results. 
The quark has dressed mass and four- momentum pi . 
The diquark four-momentum k = P — pi is constrained 
by its on-mass-shell condition fc^ — rnj., where rag is the 
mass of the diquark. 

As discussed in Ref. I, the wave function will be pa- 
rameterized by a smooth function with no singularity at 



the quark pole p{ 



If the mass of the quark and 



diquark were greater than the mass M of the nucleon, 
there would be no pole in any case, but if vriq -\- < M 
the propagator would normally contribute a pole associ- 
ated with the possibility of free scattering of the quark 
and diquark. Since the quarks are confined, this scatter- 
ing cannot take place, and detailed calculations [18] show 
that the vertex function F automatically develops a zero 
that cancels the quark pole and gives a smooth behavior 
for the wave function. Our assumption that 4" is smooth 
regardless of the value of the quark mass is a simple way 
of including confinement without investigating the details 
of this cancellation. 

The quark-diquark interaction described by the wave 
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function is parametrized through simple scalar wave func- 
tions. The diquark can have either spin-0 (isospin-0) and 
spin-l (isospin-l) components. The isospin states of the 
quark-diquark system can be written 



V3 



1 

V6 



X 



(2) 



(3) 



where t± — i: iTy are the isospin raising and lower- 
ing operators, / = ±1/2 is the isospin of the quark (or 
nucleon) 



u (orp) x' 



d{orn), (4) 



and 



^0 = ^{ud-du) 

Co = -^{ud + duj^i, 

e ^ dd = ^{^, - t^y) . 



(5) 



The operator r • in Eq. ([3]) is to be interpreted as 
transforming the initial two-component nucleon spinor 
into a two-component quark spinor. Explicitly, for the 
proton state this operation gives 



i [2d{uu) — u{ud + du)] 



(6) 



and automatically yields the correct quark content of the 
isospin-l diquark part of the proton wave function. 

In the nucleon rest frame, we choose spin states that 
are analogs of ([2]) and ([3]). This means that the spins of 
both the quark and diquark are aligned along the (arbi- 
trarily chosen) z-axis. [In Ref. I the spin of the diquark 
was chosen to be parallel to the diquark momentum, k, 
introducing an angular dependence into the wave func- 
tion from the start. Here we construct a wave function 
with no angular dependence by defining the spin with 
respect to a fixed axis.] For a nucleon moving in an arbi- 
trary direction, the spins will be aligned along the direc- 
tion of the nucleon momentum. This means that the spin 
states must be constructed in a two-step process: (i) first 
the nucleon is transformed to momentum |P| = P along 
the z direction, and then (ii) the state is transformed 
back to the final direction P. This construction bears 
some resemblance to the construction of helicity ampli- 
tudes. A more complete discussion of how these states 
are constructed can be found in Refs. [l^, [13] . 

In this paper we will need explicit results for nucleons 
moving in the ±i direction only. The scalar diquark state 
can be written 



u(0, s) u(P, s). 



(7) 



where 0j.]^/2 spin-0 diquark state, u the Dirac 

spinor, and the arrow indicates the relativistic general- 
ization. Note that, in this case, the spin of the (isolated) 
quark is specified by the spin of the nucleon. 

To construct the spin-l diquark component of the wave 
function, we begin by considering the nucleon at rest, 
where the lower two components of its Dirac spinor are 
zero. If we choose 



(8) 



with the four components e = {st, Ca;, £2} of the three 
diquark polarization states defined as in Eq. ([5]) 



£0 = {0,0,0,1}, 



(9) 



then the upper two components of fS]) will look just like 
Q and the lower two components will be zero. The state 
in the moving frame is then obtained from ([5]) by a boost, 
and becomes 



-^75 fpu{V,s) 



(10) 



where, for boosts in the z direction, the e± polarization 
vectors are unchanged, but the longitudinal one becomes 



e^p = ^{P,0,0,Sp}, 



(11) 



where Ep -|- . Note that all the polarization 

vectors satisfy 
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AA' 



P • £A = 



(where we use the metric 500 = 1 and 



(12) 
-<5,,). The 

latter condition is the usual constraint ensuring the po- 
larization vectors have only three independent states. Fi- 
nally, we note that we chose to write pH)) in terms of e* 
instead of e to allow us to interpret pO|) as amplitude for 
an incoming nucleon and an outgoing diquark in the final 
state. 

Putting this all together, the manifestly covariant nu- 
cleon wave function is the four-component Dirac spinor 

*^(P,fc) = -i=^o(P,A:)</)?7/(P,s) 

__Lv,i(P,fc)0}75 ^^u(P,.s), (13) 

which is a sum of contributions from a spin-isospin (0,0) 
diquark and a spin-isospin (1,1) diquark and Vo,i are 
scalar functions that specify the relative shape of the two 
components. If i/'o = V'l = "01 Eq- US]) reduces to 

(14) 



in precise agreement with the symmetric nonrelativistic 
SU{2) X SU{2) wave function of the nucleon. We em- 
phasize that the combination ((/)/0s + <p](f>l) is exactly 
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symmetric under interchange of any two quarks. Note 
that, because p commutes with 75/*, 5* at satisfies the 
Dirac equation 



{M- 



x(M- /')u(P,s) = 0. (15) 



The wave functions i/; are Lorentz scalars that, by 
the Hall-Wightman theorem, can only depend on scalar 
products of their arguments, and since fc^ = and 
— are fixed, they can therefore only be a function 
oiiP-kf. We chose to express this (P — fc)^ dependence 
in terms of the dimensionless variable 



(Af - m^)2 - (P - kf 



(16) 



and take a functional form for -0 that reduces to the 
Hulthen form (difference of two Yukawa functions) in the 
nonrelativistic limit, and has an asymptotic l/Q^ depen- 
dence for large Q^, as expected from_pQCD calculations 
of the electromagnetic form factors (2 it l22j | . This form is 



Nn 



Ws(/3i+x)(/32+x) 



(17) 



where Pi, p2 are range parameters (we assume P2 > Pi) 
and A'o is a normalization constant. 

We emphasize that, in the nucleon rest frame, the wave 
function p3)) contains absolutely no angular dependence 
of any kind. In Ref. I, the diquark polarization vectors, 
which were denoted by 77, depended on the diquark mo- 
mentum k with the property r] ■ k = 0. 



III. NUCLEON ELECTROMAGNETIC FORM 
FACTORS 

The calculation of the form factors, based on Fig.[lja), 
parallels the discussion of Ref. I. The relativistic impulse 
approximation (RIA) to the nucleon current in the co- 
variant spectator theory is Hi, [H [H, [U [H, [H [Hi 



(18) 



2 la^-^g^ 



2M 



where the spectator formalism places the diquark on its 
mass-shell. Our "diquark" is actually two noninteracting 
quarks with a continuous mass distribution from 2mq to 
cxD, so what we are doing here is averaging the integral 
over this mass distribution by fixing the mass at a mean 
value equal to TOs, which becomes a parameter of the 
model. (We can do this because none of the physics de- 
pends strongly on the details of this mass distribution.) 
With this restriction the four-dimensional loop integral 



reduces to an integration over the three-momentum of 
the on-shell spectator diquark 



(Pk 



(19) 



with Eg = yjm'l+'k^ the energy of the on-shell diquark, 
and the sum is over the polarizations e of the spin-1 di- 
quarks (see below). [The scalar diquark term has no 
sum.] The RIA neglects any exchange current contri- 
butions that might be present. As a consequence of our 
definitions (jl6p and (jl7p , the momentum k may be scaled 
by the diquark mass , giving final results independent 
of ma- The factor of 3 sums up the contributions from 
the three quarks. The expression is covariant and may be 
written in any frame, but is most conveniently evaluated 
in the Breit frame with the initial (P_) and final (P+) 
four-momentum of the nucleons chosen to be 

P+ = (i?,0,0,Q/2) 

p_ = (s,o,o,-g/2) 

g= (0,0,0, Q), 

the transferred four-momentum and 



= a/ -I- The spin indices of the nucleons 

have been suppressed. 

The electromagnetic coupling of a spin 1/2 quark with 
a photon is written 



3i 



- 3i\l 



g2 



J2- 



2M 



(20) 



where the subtraction term proportional to ^ is zero 
for the elastic form factors, but ensures that the current is 
automatically conserved in the deep inelastic limit. The 
use of this term for DIS was justified in Ref. ^2^, and will 
be discussed in more detail in future work. The functions 
ji and j2 are operators in the quark isospin space. For 
i = 1,2, 



(21) 



where fi± are the isoscalar and isovector combinations, 
related to the u and d quark form factors by 



If - if 

■ifid = hfi- 



^Jia — 6'''+ 2^^~ ■ (22) 

The form factors are normalized (with a — {u, d\) to 



/la(0) = 1 

/i±(0) = l 



/2a (0) = Ka 
/2±(0) = «± 



(23) 



where nu and are the u and d quark anomalous mag- 
netic moments (scaled by the quark charges) and 



^4- — t^d 



(24) 
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Note that k+ = k_ imphcs that = Hd- In aU models 
presented here, we construct fi± and /2± to satisfy the 
conditions 



hm /i±(Q2) ^ A>1 
lim /2±(Q2) ^ 0. 



(25) 



The role of A will be discussed later. The quark form fac- 
tors (j22p parameterize the charge and magnetic structure 
of the u and d CQ. 

The integral is evaluated by substituting the nu- 
cleon states (IT^ and summing over diquark polariza- 
tions. Since the diquark is a free particle, its polariza- 
tion state cannot be changed by the interaction with the 
quark, so there is no coupling between the scalar and 
vector diquarks, giving an expression of the form 



(26) 



J3T 



2M 



2M 



where -ij)^^ = ipo,iiP±, k) (and we allow for the possibility 
that i/jQ V'l for the time being) and the isospin sum has 
been done, giving j(i+2) = ^TjjiTj = - ira/^-, for 
i = 1,2. The operator D"'^ results from the sum of the 
vector polarizations of the diquarks. 

The derivation of the operator D'^'^ requires careful 
discussion. In the first version of this paper [23| our 
derivation of this polarization sum was carried out in 
the Breit frame and Kvinikhidze and Miller [28| claimed 
it was not covariant. A new, more complete and careful 
derivation of D^'' is given in Refs. [HQ, where we show 
in detail that D^"^ is both covariant and uniquely defined. 
The derivation involves the introduction of new diquark 
(and quark) polarization states which can be described 
as "covariant fixed axis polarization states." Discussion 
of these states requires the kind of care used in the origi- 
nal derivation of helicity states given by Jacob and Wick 
[29| . A general discussion of the use of covariant fixed 
axis polarization states is planned for future work. 

Using these states, the covariant polarization sum D^'^ 

is 



AP+=AP_ 
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p^P'^ ^ {p+ + p.Y{p+ + p^] 

M2 +^ 



4M2 + Q2 



(27) 



Note that this function I?^'' has no angular dependence, 
and the nuclear current is pure S-wave when — 0. 
At nonzero an angular dependence emerges from the 
wave functions , but this is due to the distortion under 
the boost and is not associated with the intrinsic struc- 
ture of the state. 

The two form factors Fi and F2 can be separated from 
the expression (|26|) using (|27|) . We now impose the con- 
dition ipQ — — -ijj as required by the symmetry of the 



state (this simple condition replaces Eq. (10) needed in 
Ref. I). Instead of reporting Fi and F2, we give the charge 
and magnetic combinations: 

Ge{Q^)^Fi{Q^)-tF2{Q^) 

= 5^(g'){ (/i+ + ra/i-) - T (/2+ + T3/2-) } (28) 
Gm{Q^)^Fi{Q^)+F2{Q^) 

= ii?(g2){ (/i+ + 5r3/i_) + (/2+ + 5T3/2-) } , (29) 

where r = Q'^/{AM'^) is not to be confused with t^. In 
these expressions the dependence of the quark form 
factors has been suppressed, and we have introduced the 
covariant nuclcon body form factor 



(30) 



The simple factorization into a product of a body form 
factor times combinations of quark form factors is pos- 
sible because only the scalar wave functions depend on 
the integration variable; the rest of the integrand (the ji 
functions and their coefiicients) are functions of only. 
We emphasize that this calculation of the form factors is 
manifestly covariant at every step. 

At Q2 = the charge form factor reduces to 

Gi=;(0) = i[l + r3]B(0) = i[l + r3] , (31) 
if we impose the normalization condition 



3(0) = [ 



ij/ = 1. 



(32) 



This equality fixes Ni^. Our theory is not complete 
enough to automatically fix the normalization (as in 
Refs. m, HI]), but we do not regard iVp as an addi- 
tional parameter. Note that the conditions ((25)) ensure 
that both Ge and Gm have the correct asymptotic be- 
havior ( ~ l/Q^ times logarithm corrections), provided 
B{Q^) ^ l/Q"* (times logarithm corrections), and this 
is easily achieved through a judicious choice of the high 
momentum dependence of the scalar wave function -0. 



IV. FIRST PREDICTIONS FOR THE FORM 
FACTORS 

Some of the predictions of this model are easily ob- 
tained from the formulas and (|29p without fitting or 
detailed analysis. 



A. Magnetic moments 

The magnetic moments are given by 

^JLp = 1 + \[k+ + 5k_) 
Mn = -| + ^(k+ - 5k_) 



(33) 
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The simplest assumption that = k_ = k gives imme- 
diately the well known quark model relation 



3 

2' 



(34) 



but the proton moment (for example) will not be even 
approximately correct unless k ^ 2 



1 



(35) 



This implies that the u and d anomalous moments are 
also equal to 2. A physical argument for this result will 
be given in Sec. IIVD] below. 

In this work we choose to reproduce the nucleon mag- 
netic moments exactly, giving 



3(/ip -f /!„) - 



1.639 
1.823, 



Ku = 1.778 Kd 1.915 



Neutron form factor 



(36) 



(37) 



Equation ([^5]) shows that Geu is identically zero if 
/i+ = /i_ and /2+ = /2- (a consequence of vector dom- 
inance if they could be approximated by a single pole at 
either the p or lu mass). Since k+ 7^ k_, Geu is not 
identically zero, but is very small. 

There are three possible ways to increase Geu- One 
possibility is to add angular momentum components to 
the nucleon. This will be discussed in a subsequent paper. 
Two remaining possibilities are (i) to break the fi+ = 
/i_ equality, or (ii) to add a pion cloud term. These will 
be discussed in Sec. IVll 



C. No scaling for Gsp and Gj\/p 

The ratio of Gsp and G]\ip predicted by Eqs. (pS)) and 
(I29p is a simple result independent of the body form fac- 
tor and dependent only on the quark form factors 



R 



G 



(/i+ + /i-)-r(/2+ + /2-) 

Gmp "(/i++5/i_) + (/2++5/2-) 



Ep 



(38) 



This will not give scaling unless the /2's are much smaller 
that the /I's, but the relations ([36)1 and ([37)1 ensure that 
this cannot be true near = 0. In fact, Eq. ([55]) shows 
that a strong violation of scaling is the most "natural" 
result. In the naive case where /i+ — /i_ = /i and 
/2+ = /2- = h = nfi, R simphfies to 



R 



G 



Ep 



G 



Mp 



fl - T/2 
/1 + /2 



1 



1 



(39) 



predicting that R will be zero at ~ 2 (GeV)^ if k ~ 2. 
However, the correct asymptotic behavior for Ge and 




FIG. 2: (Color on line) Feynman diagram for for the lowest order 
contribution to the CQ anomalous moment. 



Gm cannot be reproduced unless the /2's go to zero faster 
at high than the /I's, so the ratio R must flatten out 
at high Q^. Still, it is not hard to understand why this 
simple model predicts the violation of scaling. 



D. Angular momentum theorem 

In this section we will use the term angular momen- 
tum theorem (AMT) to refer to the light-cone result pre- 
dicting that the Pauli form factor i^2 = if there are 
no angular momentum components in the wave function 
(L = 0). 

The derivation of the AMT is based on the view of 
the nucleon illustrated in Fig. [Ijb). In this view the 
contribution shown in the figure is the matrix element 
of the bare quark current between the Fock component 
with 3 quarks and n gluons (only n = 1 is shown in the 
figure), which we represent by 

f:{q) = {P+\3q + ng) {q\r\q) (3q + ng\P^) . (40) 

Alternatively, the contributions to the nucleon form fac- 
tor can be organized as shown in Fig.[lja), where this di- 
agram is part of the matrix element of the dressed quark 
current (the CQ current) between the valence component 
of the wave function. The equivalence of these two ideas 
is represented by the identity 



light-front 



Tl = 



(41) 



(P+|3(7) (3q|P_) spectator 



with j'j the dressed quark current introduced in Eq. (|20p . 
In what follows we simplify the discussion by focusing on 
the anomalous moment of the quarks, related to the value 
of the Pauli form factors of at = 0. 

If the anomalous moments of the quarks are zero, 
Eq. ([38]) predicts that R — 1, implying that F2 — Q- 
The nonzero value of F2 in the 5'-state model comes di- 
rectly from the nonzero anomalous moment of the quark, 
and since the light-cone formalism allows for no quark 
anomalous moment, there seems to be a contradiction. 

This contradiction is easily removed by considering 
Fig. [T] as one of the many terms contributing to Eq. (|4ip . 
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In light-front language, this term is one contribution to 
the Fock matrix element ((?) involving the bare quark 
current and the Fock state with one gluon. In the specta- 
tor language, this term contributes to the matrix element 
of the lowest order anomalous magnetic moment, illus- 
trated in Fig. [21 From the spectator viewpoint, this term 
contributes to the electromagnetic structure of the con- 
stituent quark and has nothing to do with the structure 
of the nucleon itself. 

How exactly does the anomalous moment arise from 
Fig. [21)? Since the anomalous moment is the coefficient 
in the u^'^qu term, which is linear in it cannot be com- 
puted by evaluating the diagram at q = 0. Furthermore, 
a term linear in q can emerge only from the numerator of 
the diagram. Part of this numerator includes the factors 

= j d^k {mq+ i>^~ }i) 7^ {mq+ ^- /(f) (42) 

where p± = p zt ^q. As it turns out, the term linear 
in q comes from the momentum-dependent cross terms 
involving the operators /t, and hence require the 

angular- dependent terms included in a quark propaga- 
tor. If this contribution is regarded as part of the wave 
function of the nucleon, then L > components are in- 
deed required to generate a nonzero F2 , and the AMT is 
proved. However, in the spectator formalism where the 
structure of the quark current is factored out of the wave 
function, L > components are not required. 

How large is the anomalous moment predicted by the 
diagram in Fig. [2p A precise calculation would require 
summation of all the QCD diagrams taking into account 
the fact that the incoming and outgoing quarks are off- 
shell. To get a rough estimate, use the well know result 
from QED, assume that the average {p±) = m^, and 
multiply by a factor of three to include color. This gives 
(remembering that we have defined k in nuclear magne- 
tons) 



If we assume as — 1 and M/rUq ~ 3, we obtain the 
estimate k ~ 1.5. This estimate is admittedly crude, 
but perhaps sufficient to show that our phenomenological 
results for k are not unreasonable. 

What does it matter if anomalous moment contribu- 
tions are considered part of the wave function (as they 
are in the light-cone) or part of a CQ form factor (as 
they are in the spectator theory)? As long as one recog- 
nizes where various effects are included in a calculation, 
it may not. Unfortunately, our intuition is often shaped 
by hidden assumptions (which are usually based on our 
understanding of nonrelativistic quantum mechanics) so 
that perfectly ordinary results stated in one formalism 
may seem surprising when looked at from another. An- 
gular momentum components are one of these - success- 
ful nonrelativistic (or semirelativistic) models based on 
pure S'-state wave functions have existed for a long time, 




FIG. 3: (Color on line) Feynman diagram for DIS. All of the 
intermediate quarks are on shell. 

and it is important to realize that these are not "wrong" 
because they contain no angular momentum components 
in the wave function. One of the objectives of this pa- 
per is to show that these nonrelativistic models can be 
generalized to a manifestly covariant form without losing 
the intuition gained from nonrelativistic physics. 

E. Where are the lower components of the wave 
function? 

An essential feature of this model is that the nucleon 
wave function, which is the product of the quark propa- 
gator and the nucleon vertex function, Eq. (|T1), is mod- 
eled by the sum of two simple, largely spin-independent 
scalar wave functions, as written in Eq. (|13p . As men- 
tioned already in Sec. [Tll removing the quark pole from 
consideration is a simple way to build quark confinement 
into the model, but removing the spin-dependent posi- 
tive energy quark projection operator from the nucleon 
wave function has the effect of removing the lower compo- 
nents of the quark spinor and replacing them by the lower 
components of the nucleon spinor. Since these lower com- 
ponents are exactly zero when the nucleon is at rest, the 
bound quark must also have no lower components, and it 
is natural to ask what kind of interaction could produce 
such a result. 

To answer this question, consider the covariant specta- 
tor equation that might describe the nucleon as a bound 
state of an off-shell quark and on on-shell "diquark" 

S-'iP~p)^{p,P)^ f V{p,k;P)^ik,P) (44) 

Jk 

where J^, was defined in ^ is the wave function 

defined in Eq. (|T1), and S^^ is the inverse of the quark 
propagator 

S~\pi) ^ mg- . (45) 

A kernel V can be constructed that will give linear con- 
finement in coordinate space [l^, [s^l and the details of 
its construction need not be discussed here. In the non- 
relativistic limit (defined by nig 00) Eq. (|44| has the 
nice feature that it reduces to the Dirac equation 

(i ^-mq-Vir)j^{r) = 0, (46) 

where \E'(r) and V{r) are three-dimensional Fourier trans- 
forms of their momentum space counterparts. 
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To obtain a wave function without lower components 
(in the rest frame) it is sufficient to replace the kernel by 

Vip, k- P) = 0(p, P)V{p, k- P)0{k, P) (47) 

where 0{p, P) = S^^{pi)[M+ pi = P—p, and for any 
Dirac operator A = ^^A^j'^. The construction preserves 
the hermiticity property of the kernel, V — V. With this 
substitution the wave function satisfies a new spectator 
equation 

^{p, P) = [M+ p] [ V{p, k; P)0{k, P)^{k, P) (48) 

Jk 

showing immediately that 5* satisfies the Dirac equation 
psp . sufficient to ensure that its lower components are 
zero in the nucleon rest frame. 

This demonstration shows precisely how the features 
of the spectator theory can be exploited to build in this 
phenomenology: since the quark (in this case) is off-mass- 
shell, and since the kernel V{p, k; P) always depends on 



the total four-momentum, the substitution (j47p does not 
violate any of the constraints that must be satisfied in 
constructing a phenomenological kernel. 



V. DEEP INELASTIC SCATTERING 

Before studying the details of the predictions for the 
form factors, we will show how this model is able to de- 
scribe deep inelastic scattering (DIS). 



A. Quark distribution function 

The DIS cross section can be calculated from the imag- 
inary part of the forward handbag diagram shown in 
Fig. [21 The cross section depends on the hadronic tensor 



. -4,./ { (»: - if) - (p„ - ^) [p. - ^) («) 



where Cq — y "^g + P'^ the energy of the on-shell quark 
in the final state. Eg is the diquark energy encountered 
before, Wi and W2 are the DIS structure functions, and 
is the hadronic current 



= -u(p',Si)j^(p',p)*Ar(P,fc) 



(50) 



with Jiiip'^p) the quark current defined in Eq. (|20p and 
^jv the nucleon wave function defined in Eq. p^ . Here 
the isospin / and the polarizations of the diquark and 
nucleon, §2 and s, have been suppressed. 

Detailed evaluation of the hadronic tensor is sketched 
in Appendix A. After the spin sums have been carried 
out, the structure functions become 

Wi = ^ ff {P-p')^'iP,k) 



2ttM 



p' k 



= ^ // 2M^Aij^iP,k), (51) 



p' k 



2ttM 



where A was defined in Eq. ([25|) . the isospin-dependent 
charge operator is 



■T3 



2et 



2e? 



proton 
neutron 



(52) 



the integral is 



Ilp'k II 



d^p'd^k , 



(2^)2 AeqEs 
d^k 



6^ip' + k~P~q) 



5+{m'-p")5+{mi-k'), (53) 



and 



A = 



(27r)2 "^^ 1 

jp' ■ P)Q^ + {P ■ q)ip' ■ q) 
NPQ^ + (P • g)2 



(54) 



The expressions for Wi and W2 are covariant, but it is 
convenient to evaluate them in the laboratory system us- 
ing light-cone coordinates. In our notation, an arbitrary 
four-vector v is written 

w = v^} v±=vo±vs, (55) 

so that the scalar product is 

V ■ u — v^u^ = i(ti+u_ + t^-U+) — • u_L (56) 

The four vectors in the laboratory frame are 



P = {M, 0,0,0} 



q 



2Mx 



0,0, x Q 



(57) 
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where x is the usual Bjorken scaUng variable. In the DIS 
limit, in light-cone form, these momenta become 



P = {M,Af,0} 



(58) 



Using this notation, we find that the structure func- 
tions Wi and W2 scale, and satisfy the Callen-Gross re- 
lation 

VW2 = 2MxWi EE e] xfg{x) . (59) 
with a quark distribution function fq{x) given by 



~ Att J (27r)2(l 



-V'^(XDIS) 



(60) 



where xdis is the value of x, Eq- ifTB]). in the DIS hmit. 
If 



and k\ — Mrrisy, then 



XDIS = -2 



M 



r + y ^ (1 - x) 
I — X r 



(61) 



(62) 



The values of k± were fixed by the S function restraints in 
Eq. (|53p . leaving an integral only over the perpendicular 
components of k. Using the wave function (|17p the k± 
integral is easily evaluated, giving 



where 



H = 



16^V(/3i - (32)' 



(63) 



1 



1 



2 R2 





r 1 - 


A-2- 


f + 

1 — X r 


X 




Mrris 


x 1 — a; 



(64) 



and = Mms/3j + (M — ms)^. 

Before we examine the implication of Eq. (|63p . we 
study the normalization of fq{x). 



Multiplying both sides of (|65p by A^, defining fc_ = 
M(l — x), and using ((63)) gives the relation 

= / dx fq{x) (spectator theory) . (66) 
./ — 00 

The contributions to this normalization integral from 
the region x < are a feature of the spectator theory and 
arise because, by design, the spectator theory includes 
only contributions from the singularities of the spectator 
[Slj . If additional singularities in the wave function are 
also taken into account, the contribution to the norm (j65p 
from the region of x between —00 and will be canceled. 
In Appendix B we show that including the singularities 
of the wave function limits the fc_ integral to the region 
[0,Af], giving 



A^ = 



dx fq (x) (light- front) . 



(67) 



Hence both theories (spectator and light-front) give the 
same quark distribution amplitude; the only difference is 
the normalization condition ([66]) or (|67p . 

To compare with experiment, we use the valence dis- 
tributions obtained from a recent global fit [l^l- Results 
will be compared to the proton average 



cxpt 



4 

gXUV 



hxd\ 



(68) 



(where xuy and xdy are Eqs. (1) and (2) of Ref. [35 
and is normalized to 



dxfqix)U^pt = 1.006 - 1.00 



dxxfq{x)\cKpt = 0.171 



(69) 



For the light-front, consistency therefore requires choos- 
ing A = 1, while for the spectator we must have A > 1 
so the integral over the physical range < x < 1 will be 
exactly unity. 

We now turn to a discussion of the fits to the form 
factors and the quark distribution functions. 



VI. RESULTS 



B. Normalization of the quark distribution 
function 



The normalization condition ([32]) for the nucleon wave 
function can be expressed as a normalization condition 
for the quark distribution function. Writing the normal- 
ization integral in terms of light-cone variables 



1 = 



dk^ dk^ d'^k± 

r°° dk^ r d^k 

^- 



2(27r 



6+{ml-k^)^^P,k) 



(65) 



A. Vector dominance model for the quark form 
factors 



Motivated by the vector dominance model (VDM) , the 
quark form factors are written in the following form 



/i±(Q') - A + 



(1-A) 



c±Ql/Mi 



1 + Ql/ml (l + Ql/Mly 
d± , (1 - d±) 



1 + Ql/rnl 1 + Ql/Ml j ' 



(70) 



where is the p (w) mass for the isovector (isoscalar) 
form factors, and = 2 is a fixed heavy mass intended 
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FIG. 4: (Color on line) The first few terms of the bubble sum for the quark form factors. The first term is the bare charge, missing from 
the /2 form factors. 



to approximate the contributions from all heavy vector 
mesons in each channel (this could depend on the isospin 
but this possible dependence is ignored in these fits) . All 
masses and momenta are in units of the nucleon mass M, 
so that, for example, Qq — /AP . 

It may be helpful to give a brief justification for the 
quark form factor models ([70)1 . From a microscopic point 
of view, the VDM for the form factors emerges from a 
sum similar to that shown in Fig. 3) Here, the quark- 
quark interaction is approximated by a contact term, giv- 
ing 



Model j3i,(52 c+,C- d+,d- 6i5,&j\/ A, r Nq,x^ 



h = 



Xe + gB{s)\e ^ 
g\eB{s) 



[gB{s)Y\e 



Xe -I- 



l^gB{s) 



(71) 



Choosing a suitable B{s) that goes like Q ^ as 



oo then gives Eq. ([70)1 for the form factors fi±. The 
construction of /2 is similar, with the first term missing 
because there is no bare anomalous moment at infinite 



1(4) 0.057 


2.06 


-0.444 




1.22 


10.87 


0.654 


2.06* 


-0.444* 




0.88 


9.26 


11(5) 0.049 


4.16 


-0.686 




1.21 


11.27 


0.717 


1.16 


-0.686* 




0.87 


1.36 


111(6) 0.078 


1.91 


-0.319 


0.163 


1.27 


12.36 


0.598 


1.91* 


-0.319* 


0.311 


0.89 


1.85 


IV(9) 0.086 


4.48 


-0.134 


0.079 


1.25 


8.46 


0.443 


2.45 


-0.513 


0.259 


0.89 


1.03 



TABLE I: The two lines for every model give the values of the 
8 possible (9 for Model IV) adjustable parameters and the 4 
constants fixed by the constraints. The fixed constants are A, 
the asymptotic value of /i given in Eq. (|25p . the diquark mass 
ratio r, Eq. (|61|l . the normalization parameter Aq, and the 
per data point. In each case k+ = 1.639 and K- = 1.823. 
Parameters labeled with an * were constrained during the fit 
to equal the one above it. For Model IV the heavy mass 
defined in Eq. (|70|) was also adjusted, giving a best value of 
Mh = 2.556. 



B. Fitting procedure 

During each step in the fit, we first normalized the 
wave function using Eq. (I32|) . determined the anomalous 
moments of the quarks by fixing the proton and neutron 
magnetic moments, and then determined the asymptotic 
charge parameter A and diquark mass vris by satisfying 
the following two constraints 

(/) ^ Cdxf,{x) = l 
Jo 

(xf) = [ dxxf.ix)^ 0.171 (72) 
Jo 

where the quark distribution amplitude fq{x) was defined 
in Eq. (|60|) . Even though the form factors are indepen- 
dent of fUs (recall that it was scaled out of the results) 
the deep inelastic quantities are not. Since fq includes A 
in its definition, to^ was first determined from the ratio 
{xf) I (/) , and then A was fixed by the normalization re- 
quirement (/) = 1. This procedure was quite stable and 
led quickly to a good determination of the other param- 
eters. 

In this paper we chose to define the charge of the u 
quark (for example) to be 2/3 at Q-^ = 0, and to nor- 
malize the wave function to unity. QCD suggests that it 
might be preferable to define the u quark charge to be 
2/3 at = oo. This can be done by dividing the quark 



I II III IV expt 

rl 0.764 0.791 0.851 0.703 0.780(25) 

rl -0.005 -0.104 -0.241 -0.102 -0.113(7) 

X?adii 118 0.89 170 6.11 



TABLE II; The values of the charge radii and the xLdii foi^ each 
model. 

form factors by A which changes the charge at = to 
2/(3A). If the wave function is then normalized to A (in- 
stead of unity, as is now the case), all of the formula for 
the form factors will be unchanged. However, the numer- 
ical results will change slightly because the normalization 
requirement ([7^ will change to (/) = A (keeping the def- 
inition of ]q unchanged). This will give the same value 
of 771s, but a slightly different value of A, leading to a dif- 
ferent minimum with slightly different numerical results. 
We save discussion of this subtlety for future work. 



C. Fits without a two-pion cut 

Since the VDM for the form factors includes a p pole 
and the p couples to the 27r channel, the quark form 
factor includes two-pion terms. It is therefore next to 
impossible to separate out a pion cloud term from the 
quark form factor contributions. However, the two-pion 
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FIG. 5: (Color on line) Data for the nucleon form factors compared with the four models discussed in the text: Model I (dotted line), II 
(short dashed line), III (long dashed line) and IV (solid line). The Gjv/p data is from Arrington | 3>^ . and the Gjv/n data includes that used 
by Bosted (solid circles) in his global fits to the form factors [3j[ and the unpublished data from JLab Hall B (open circles) |35j . The Gep 
data are from JLab Hall A, Jones ei. aZ.,and Punjabi et. ah, [l|, \§i (squares) and Gayou, et. ah, Ql (triangles). The Geu data are single 
Q2 points from MAMI [H, [H, NIKHEF [H, and MIT-Bates [3?| (solid circles), and from JLab Hall C by Zhu, et. al.,]^ (triangle). 
Warren, et. ah, [4l| (squares), and Madey, et. al., j42j (diamonds), and from MAMI by Glazier [J^Wopen circles). Only Gev. data obtained 
from deuteron targets are included. For a list of the data see the nucleon form factor data base |44| . 



channel contributes a cut starting at = — > 4m J, 
and it has been known for many years that the singular- 
ity associated with the onset of this cut, which goes like 
(g^ — 4to^)^/^, has a strong effect on the form factors near 
— 0. The same cannot be said for the isoscalar form 
factors. The isoscalar singularity is produced by the on- 
set of a three-pion continuum, with a three-body phase 
space starting smoothly at 9to^. It is much weaker, more 
distant, and less important. 

The VDM model described in Eq. (|7D|) does not have 
the two-pion singularity. In the next section we will dis- 
cuss how this singularity can be added, but here we look 
at the consequences of ignoring this singularity. 



We present four models with different features that 
represent different aspects of the physics. The parame- 
ters for each model are shown in Table HI In this sec- 
tion we discuss the first two of these models, I and II. 
These models do not contain the pion-cut singularity. 
Model I is simplest model with the fewest parameters 
(4). It assumes approximate isospin symmetry by con- 
straining /i+ ~ /i_ and /2+ ~ /2- (they are not equal 
because nip ^ and k+ 7^ As already discussed 

in Sec. IIVBI approximate isospin symmetry for a pure 
S'-wave model of the nucleon gives a very small neutron 
electric form factor, and this is show in Fig. [5] and Table 
im (showing that the charge radius of the neutron is off by 
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FIG. 6: (Color on line) Log plot of the same curves and data shown in Fig. [5] [For an explanation sec the caption to Fig. [5]] These curves 
show the low structure discussed in Ref. [45j . 



almost two orders of magnitude) . We do not believe that 
Model I gives a credible description of the form factors, 
even though its predictions are very close to many other 
simple quark model calculations. We present it to show 
what can be achieved from very few assumptions. 

As discussed above, one of the ways to obtain a reason- 
able neutron charge form factor is to break the isospin 
symmetry. Model II allows the isospin to be maximally 
broken by the /i form factors, but preserves the approxi- 
mate isospin symmetry of the ji form factors. This adds 
only one more parameter, but does an excellent job of 
fitting the data, and incidentally, gives almost perfect 
results for the two charge radii. This model is perhaps 
the best compromise between the need to obtain a good 
description of the form factors, and the phenomenologi- 
cal requirement that the description be economical, using 
only a few parameters. 



Now we discuss the effect of adding a term that explic- 
itly includes the two-pion threshold singularity. 



D. Including a pion cloud (two-pion cut) 

The imaginary part of the isovector form factors has 
been calculated in chiral perturbation theory (ChPT) [i^, 
firj . and also extracted from data [1^. We model this 
behavior by adding extra terms to Ge- and Gu- of the 
form 



AG,_(Q2) = 



6,; 1 



r(Qf) 



(73) 
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FIG. 7: (Color on line) Im GE—ft^ compared to the experimental 
fit of Ref. \4m (upper solid curve) and the one loop ChPT calcula- 
tion of Ref. [46l | (lower dotted curve ending at t = 20). The units 
are pion masses, p. = m-,^. The dashed lines, in order of decreas- 
ing size, are the models with values of feg = 0.55, 0.45, 0.30, and 
0.20. The two dot dashed lines correspond to Model III (upper) 
and Model IV (lower). 



where the width function is 



(74) 



The choice of this function is discussed in some detail in 
Appendix C. The parameters a^, mi and TO2 were cho- 
sen to reproduce the threshold behavior of the imaginary 
parts of the isovector Ge and Gm form factors, calcu- 
lated in chiral perturbation theory, and, for Ge, also ob- 
tained from experiment. The imaginary part near thresh- 
old is insensitive to the value of the additional parameter 
hi adjusted during the fit. The behavior of the imaginary 
parts is shown in Figs. [7] and [S] 

Models III and IV both include the pion cloud term 
([73)1 . Model III shows that adding this term signifi- 
cantly improves the description even when the approx- 
imate isospin symmetry used in Model I is still main- 
tained. We conclude that we may significantly improve 
the fit to the data by either breaking the symmetry of 
/i or by adding a pion cloud term. However, these two 
alternatives are not completely equivalent. Breaking the 
symmetry of /i (Model II) produces a better fit (lower 
X^), much better radii, and also gives a Geu with a bet- 
ter shape. The shape of Geu seems to be particularly 
sensitive to how the isospin symmetry is broken. 

Finally, Model IV shows that a combination of the two 
effects studied separately in Models II and III can give 
a precision fit to the data. [To obtain this the heavy 
mass Mh was also varied, adding one more parameter.] 
The expected for a statistically perfect fit should lie 



in the range x = 1 =t 



where n is the number 




of data. Here the number of data is 117, leading to an 
expected x^/n = 1±0.13 [i^, nicely consistent with our 



FIG. 8: (Color on line) Im Gm-I^^ compared to one loop ChPT 
calculation of Ref. O (lower dotted curve ending at f = 20) . The 
dashed lines, in order of decreasing size, are the models with values 
of bjv/ = 1.50, 1.00, and 0.55. The two dot dashed lines correspond 
to Model III (upper) and Model IV (lower). 



result of 1.03, or 1.12 if the radii are included in the fit. 
One of the reasons a high quality fit is possible is that we 
used data for Gmp recently reanalyzed by Arrington [3^ ; 
these data are more consistent with the new Ge^ data, 
and this consistency is required for a good fit. The new 
JLab high precision Gun data are not yet final [35] and 
were therefore not included in our fit; if the final data set 
remains close to the preliminary one a refit of the model 
may be necessary, and we have not investigated the effect 
it will have. The model predicts that Ge^ vanishes for 
Q2 0.7.5 GeV^. 



E. Quark distribution functions 

The results for the quark distribution function ob- 
tained from the four models are shown in Fig. [9] Note 
that all the models give very similar results; the distribu- 
tions are all too narrow and approach zero as a; ^ 1 too 
fast, but at least give a qualitatively reasonable descrip- 
tion. Better agreement would have been obtained if we 
had compared the models with a quark distribution func- 
tion evolved to a higher (4 GeV^, for example), which 
would be somewhat more peaked than the one shown in 
the figure (evolved Xo = \ GeV'^). We have made 
no attempt here to improve the agreement between our 
result and experimentally determined distributions; this 
is a subject for future work. 

Finally, Fig. 1101 shows the normalized momentum dis- 
tribution for the nucleon in its rest frame (where it is 
spherical) . The larger size of the wave functions for Mod- 
els I and II are compensated by a slightly smaller size at 
larger fc; the difference is exaggerated in the figure be- 
cause the factor A;^, included in the normalization inte- 
gral, is excluded from the distribution. 
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FIG. 9: The valence quark distribution xfq{x) as a function of 
X. Each curve is normalized to 0.171, the experimental momentum 
fraction carried by valence quarks in the pro ton. The thick short 
dashed line is the empirical fit from Ref. [33 | (evolved to = 1 
Ge V'^ ) , and the other lines (with the same line style used in Fig. [5] 
and nearly indistinguishable from one another) are the four models 
discussed in the text. 
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FIG. 10: The distribution p(A;) = mf|V'o(P,fc)|^/(2£;s) as a func- 
tion of the momentum k for each of the four models discussed in 
the text (drawn with the same line style used in Fig. |5j. This 
distribution is normalized to An Jg°° k'^dk p{k)/{2iTm,s)'^ = 1. 



VII. SUMMARY AND CONCLUSIONS 
A. Overview and comparison with previous work 

This work presents a simple covariant model for the 
nucleon based on the foUowing assumptions: (i) the nu- 
cleon is composed of three valence constituent quarks 
(CQ, massive extended particles dressed by the quark- 
antiquark interaction, the pion cloud, and gluon sea, all 
parametrized by quark form factors), (ii) the three-quark 
system is described by an internal wave function consis- 
tent with the properly symmetrized covariant spectator 
formalism and the overall center-of-mass motion of 
the total system is described by a free Dirac equation for 
a particle of mass M (the nucleon mass), (iii) the internal 
wave function has a structure built entirely from 5-wave 



components with exactly the same spin-isospin content 
as the simplest nonrelativistic SU{2) x SU{2) model, 
(iv) the CQ form factors are normalized to the quark 
charges in the confinement limit (Q^ = 0) and reproduce 
the behavior of pointlike quarks in the large regime 
giving form factors that behave (up to logarithms) like 
pQCD at very large Q^, and (v) the wave function has 
the appropriate nonrelativistic limit. A nice feature of 
this model is that the two components of the wave func- 
tion, corresponding to spin-0 and spin-1 diquark states, 
are described by one scalar wave function, as required by 
symmetry under quark interchange. 

Aspects of this model have been included in previous 
work. The use of CQ form factors is certainly not new, 
dating back to the early valon mod el ISOll. and has been 
pursued in many more recent papers [5ll |52| | . Others have 
used a vector dominance model to describe the form fac- 
tors, with an additional phenomenological factor added 
to simulate the body form factor (ssl . Im . Issj . The two- 
pion cut is also included in this work [53| . Following dif- 
ferent lines, pQCD inspired calculations have appeared 
m, US] , as well as other calculations based on general- 
ized parton distributions [111, QCD sum rules [5§| and 
lattice QCD [60|. A comprehensive review can be found 
in Refs. [7l, [Sj, [61|] . 

The fits presented in this paper were obtained by first 
adjusting the quark anomalous moments so that the mag- 
netic moments of both the neutron and proton are repro- 
duced exactly, and then adjusting the "diquark" mass 
and the asymptotic quark charge to give the correct ex- 
perimental result for the number and total momentum 
fraction carried by valence quarks in DIS. With aspects 
of both the low and high energy behavior fixed, it is then 
possible to study the sensitivity of the form factors to 
various physical assumptions. 

The physics is illustrated using four models of increas- 
ing complexity. Model 1 examines the consequences of the 
simplest assumption: the quark form factors are approx- 
imately independent of isospin (some breaking is built in 
because mp ^ and k+ 7^ k_). With only four param- 
eters (two in the wave function and one each for /i and 
72) the dotted lines shown in Figs. [H] and [H] are obtained, 
giving a good description of all the form factors except 
for the neutron charge form factor, Geti- Similar results 
were obtained in recent work by Julia-Di'az, Riska, and 
Coester [g^l- Note that the nonscaling behavior of the 
proton charge form factor is easily described. Ironically, 
Geu is far more important to our understanding than 
the more popular Gep measurements. 

To describe the neutron charge form factor in a pure 
S-wave model it is necessary to break the isospin invari- 
ance. There are two ways to do this: split at least one 
of the quark form factors (/i+ 7^ /i_ or /2+ 7^ /2-), 
or add a term describing the two-pion cut (starting at 
t = = —Q^ > 4toJ). These two methods are, to 
some extent, merely different ways to model the same 
physics: the two-pion cut coming from the p contribution 
(isospin^l) is much stronger and has more structure at 
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low t than the corresponding three-pion cut associated 
with the uj contribution (isospin=0). This mechanism 
automatically breaks isospin symmetry. This feature has 
been included in the recent vector dominance models of 
Bijker et. al., \5^ . 

We found that breaking the /i symmetry was far more 
efficient than breaking the /2 symmetry, and Model II 
shows the dramatic effect of allowing c+ and c_ [defined 
in Eq. ([70)1 ] to vary independently. 

Model III achieves a similar result, but its charge radii 
are quite poor, and the Geu given by this model has the 
wrong shape. 

Finally, Model IV shows that we can obtain a precision 
fit if both the symmetry of fi and /2 is broken, and the 
pion-cut term is added. 

Comparing Models II, III, and IV, we conclude that we 
can fit the data by including an explicit pion cloud and/or 
adjusting the vector dominance contribution (coefficients 
c±). This result is consistent with Ref. ^63|, which states 
that "an unambiguous extraction of the pion cloud con- 
tributions is not possible." Furthermore, the pion cloud 
descriptions of Refs. tH, 113 "^o ^^t determine the real 
parts of the form factors, except possibly in the limited 
region < 0.1 GeV^ (as noted in Ref. [Sj). In any case, 
our model for the pion cloud contribution for the nucleon 
form factor in the spacelike region gives a contribution 
that differs only slightly from the charge (or magnetic 
moment) in the small region, and decreases quickly 
with increasing . We conclude that we can describe the 
nucleon form factor data using vector dominance only. 

In all the models we have very different values for the 
parameters Pi and P2 (one order of magnitude). This 
is evidence of an almost perfect separation, in dynamics, 
between the low momentum (/3i) and the high momen- 
tum {P2) regimes. 

The overall success of this family of models, each of 
which is consistent with the expected pQCD asymptotic 
behavior (up to logarithms) and describes the physics of 
DIS scattering qualitatively, shows that the form factor 
data do not necessarily demand that the nucleon wave 
function include L > nonspherical angular momentum 
components (although they are certainly not ruled out). 

In forthcoming work, including Ref. [20j, the structure 
of the nucleon wave function will be used to describe the 
transitions to other baryons and to excited states of the 
nucleon. For that purpose a modified version of the di- 
quark propagator D'^'^ will be used. These studies will 
give more definitive information about angular momen- 
tum components in both the nucleon and other baryons. 



B. Shape of the nucleon 

We now tm'n to the interesting discussion of the shape 
of the nucleon. Shapes of nuclei have been discussed for 
over 50 years, and are still an active area of research. The 
Contemporary Physics Education Project (CPEP) [gH] 
may be illustrative of the discussion found in the recent 



popular literature. In the discussion of their nuclear sci- 
ence wall chart, they refer to nuclear shape without giv- 
ing a definition. In an early paper. Rainwater 66] relates 
nuclear shape to the presence of a nuclear quadrupole 
moment, and this definition seems to be universal. A 
very nice discussion of nuclear shapes and quadrupole 
moments can be found in the Nobel lectures by Rainwa- 
ter Bohr [H, and Mottelson |[69|]. 

The quadrupole moment is determined by the charge 
density. For example, consider a charged spin- 1/2 quark 
moving about a fixed spin-1 particle (the diquark). Since 
the diquark is merely a representation of the two quarks 
not being probed by the photon (we assume that there 
are no two-body charge operators), we can ignore its 
charge. For definiteness, suppose the quark has angu- 
lar momentum £ — 1 and the total angular momentum 
of the state is 1/2. In nonrelativistic terms, the charge 
density of the spin-up state is then 

Pe(r) = e?/'^(r) 



( 5 1 ^1*1 V 3 (5 1^1 



10 



47r 



sm 



cos' 



(75) 



This is a spherical result, even though the individual com- 
ponents that make up the wave function are not spheri- 
cal. 

In our model, the momentum space charge operator at 
= (both nucleons must be at rest to avoid Lorenz 
contraction effects) can be extracted from Eq. (|28p 



Pe(k) = 3*Ar(P,A:)ji7o*jv(P,A:) 
^ i(l + r3)l^(F,fc)l^ 



(76) 



where ji , defined in Eq. ([20)1 , is the quark charge opera- 
tor. Note that Pe(k) is proportional to the nucleon charge 
operator and to the square of the scalar wave function ip 
(a function of only) as expected. The electric charge 
density of the nucleon is spherically symmetric. The ra- 
dial momentum distribution \ip(P, A;)]^, multiplied by the 
kinematic integration factor m^/{2Es), is shown in Fig. 
[TU]for the two models considered in this work. 

The charge distribution derived in Ref. I is also spheri- 
cally symmetric, even though the wave functions are not. 
This is because the angular dependence in the wave func- 
tions of Ref. I (due to the angular dependence of the po- 
larization vectors 77) is canceled in much the same way as 
in example ((75)) . 

We can also use the matter distribution to analyze the 
nucleon shape. In this case the operator ji is absent, and 



,(k) = 3^NiP,k)-fo^NiP,k), 



(77) 



This is also spherical. 
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If the shape is defined by the charge density, we con- 
clude that the nucleon is spherical, unless it can be shown 
to exhibit some collective motion that would allow us to 
interpret it as a deformed state precessine about an axis 
different from its symmetry axis [6i,[63,|6i,[6i,|73. This 
would imply the existence of a rotational band and su- 
perlarge quadrupole radiation, which have not been ob- 
served. 

Recently Miller Yl] and Kvinikhidze and Miller ^ 
introduced a spin direction dependent (SDD) density op- 
erator as a means of describing the shape of the nucleon. 
Nonrelativistically, the charge SDD operator is 



SDD 



(r,n) =pe(r)i[l + cr-n] 



(78) 



where n is the direction of the spin quantization. Choos- 
ing n to be in the +z direction, the spin-up matrix ele- 
ments of this operator for the nonrelativistic case studied 
in Eq. ^ are 



SDD/ 



n) 



( 5 1 ^1*1 V H 5 1 ^1 



5(1 + ^3) 



47r 



■ cos^ 6, 



(79) 



because, in this case, the spin projection operator has 
projected out the spin-up state of the quark, unveiling 
the angular momentum contained in Y-^q. Using this def- 
inition Kvinikhidze and Miller 13] analyzed the relativis- 
tic nucleon model of Ref. I, demonstrating that the SDD 
matrix elements generate angle- dependent terms in the 
nucleon SDD charge and matter densities. 

Relativistically, the SDD density operator includes the 
factor of (70 + n • 775)72, where n is a unit vector that 
specifies the quark spin direction. With this definition, 
the SDD electric charge distribution for a nucleon with 
polarization s becomes 



pr^(k,n,s) 



-^'Ar(P, /c)jl(70 + n • 775)1'Ar(P, fc) 



I + T3 , l + 5r3 



and for the matter distribution 



6 



-n • s 



(80) 



^SDD 



(k,n,s) = -§Ar(P,fc)(7o + n-775)*Ar(F,fc) 



2 

= i(3 + n.s)|VXF,fc)|2 



(81) 



where, as before, ipiP, k) is independent of angles in the 
nucleon rest frame. For the model presented in this pa- 
per, these distributions are also spherically symmetric, 
refiecting the fact that the model contains only S'-wave 
components. 

What are the larger implications of these observations? 
The SSD operators are potentially quite interesting. If 



they can be measured directly, they will reveal the an- 
gular momentum content of a state. But using these 
operators to define the "shape" of a state is contrary to 
what is usually understood as the shape (i.e., the charge 
or mass quadrupole moment density). 
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APPENDIX A: EVALUATION OF THE 
HADRONIC TENSOR IN THE DIS LIMIT 

Starting from the general formula (j49p . the hadronic 
tensor requires evaluation of the following integral and 
trace 



p' k 



f^^^^iP,k)[f,TTo - i(r,j\V,)Tri}(Al) 



where the DIS current was defined in (jSOp , the quark cur- 
rent ji in ([2T|) , and the traces for contributions from the 
isospin and 1 diquarks are (here we use the normaliza- 
tion uu = 2M) 



Tro = Tr^iM+n[l'-f)/[l"-f, 

Tn - iTr|(M+ nri' (7^ - ^) / 



x(7^-^j7V|i^a,, (A2) 

and the phase-space integral was defined in (|53|) . The 
quark mass has been dropped from the traces since its 
contribution is negligible in the DIS limit. This tensor 
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is covariant and may be calculated in any frame, and at 
this point we will not specify a frame. 
The quark currents are 

Jf - + ^.r,y = A^(^ + ir3) 

3^»(^ + ^^3)t. = ^-i^T3 (A3) 

The diquark polarization sum, defined for the case when 
P+ = P- +9 in Eq. is now simply -gap + PaPp /M'^ , 
from which we obtain immediately 



Tro = -Tn = 4PV^ + 4^"^^^ 
-4(P.p') (g^^^- 



where, for any four-momentum. 



(A4) 



(A5) 



Hence the quark currents (|A3p must be added, giving the 
current factor reported in Eq. (|52p. 

To extract the structure functions Wi and W2, we go to 
a collinear frame (where P and q are in the z direction) . 
Then, since the wave function depends only on (P—k)^ — 
{p' ~ 9)^1 E^nd since the numerator is linear in p', the 
transverse components of p' integrate to zero, and we 
can write 



p = AP + Bq 

with A given in Eq. ([51)1 and 

{p' ■P){P-q)-M\p' -q) 
(p ■ qY + M^Q^ 

Since q — 0, the trace terms give 
Tro = -Tri = SAP^P^" - 4(P • p') ( gf"" 



(A6) 



(A7) 



(A8) 



giving the expressions (jSip for the structure functions. 

If the integrals over k± are evaluated in the nucleon 
rest frame, as discussed in Sec. IV Al the phase-space in- 
tegral can be written 



Ilp'k 



dkj-dk- 



(1 -x)Q2 J 2(27r)2 



giving the limits 

P-p 



2x 



A 



(AlO) 



These lead directly to the expression ([60|) for the quark 
distribution amplitude. 



APPENDIX B: WAVE FUNCTION 
NORMALIZATION IN SPECTATOR AND 
LIGHT-FRONT THEORIES 

To simplify the formulas in this section, we work with 
a wave function with only one pole: 



^{P,k) 



No 



"m-siP -I- X - ie) 



(Bl) 



All results are also true if the wave function has two 
poles, as can be easily seen after the arguments have 
been developed. 

Begin by examining a four-dimensional integral related 
to the normalization integral. At = 0, the integral is 



d*k V^(Po,A:) 



(2^) 



d*k M^N- 



(27r)4 DiDl 



f (B2) 



where the spectator four-momentum k is unconstrained 
(so that the spectator is no longer on-shell) and the de- 
nominators therefore contribute two poles and two dou- 
ble poles in the complex fco plane: 



Di 
D2 



mi - k^ 
MmsP- 



+ (M - ; 
M~ko' 
IkF-M 



- ie){Es + ko — ie 
f - {M - kof - 



(B3) 



where ^ = MrusP + {M — rus)^- The single poles are 
at fco — ±Es and the double poles are at fco = M ± 

The spectator theory organizes the infinite series of 
Feynman diagrams that describe (in this case) the quark- 
diquark interaction so that the spectator (diquark) is 
always on-shell; contributions from terms in which the 
spectator is off-shell are included in higher order terms 
in the kernel. In the present example this means that the 
only contribution from the integral (jB2j that we may 
include comes from the positive energy spectator pole, 
and other contributions to the four-dimensional integral 
would be included as part of the higher order terms in the 
kernel, and hence are part of the wave function ip that we 
are modeling. Separating out the positive energy spec- 
tator pole from (jB2p gives the spectator normalization 
integral 



P 



d^k 



(27r)32P, 



■V^'(Po,A:) 



(B4) 



where fc^ = m^. 

In light-cone theory the wave function is the solution of 
a generalized hamiltonian dynamics with — H + P^ 
the generalized hamiltonian. The -I- component of mo- 
mentum is not conserved. However, some of the features 
of the light-cone approach (and, in particular, the nor- 
malization issues discussed here) can be understood by 
evaluating the integral (|B2[) using light-cone coordinates. 
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In these coordinates, the four-dimensional integral be- 



comes 



h = 



f dk+ dk^ d^kj_ A'PN§ 



(B5) 



k^k^ 



where now 

D[ = ml + h 
D'^ = £,+k\-{M -k+){M -k-)-ie. (B6) 

This integral has only one pole and one double pole in 
the complex k+ plane. If < fc_ < Af, the single pole 
is in the lower half plane and the double pole is in the 
upper half plane, but outside of this region either the 
single or the double pole migrates so that all are in the 
same half plane, and the integral is zero [UIt^I- Hence 
the exact answer is given by the spectator pole, with 
Q < k- < M. Introducing l~x — k^/M, the momentum 
fraction carried by the diquark, the integral becomes 



with 



dx 



1-xJ 2(27r)3 x^d^ 



kl 



kl 



1 



(B7) 



(B8) 



At first glance, it might appear that the light-cone in- 
tegral (|B7p has nothing to do with the spectator integral 
(|B4I) . but an interesting connections was noticed a long 
time ago. To make this connection explicit, we trans- 
form the spectator integral into light-cone variables by 
transforming kz to x using 

M{l-x) = Es~kz 

Mdx = (^l-^^dkz=M{l-x)^ 

M{Es + kz) = iml + kl)/{l-x) 

D2 = D2{k) - ^ + kl + kl-{M- Esf 

= xd. (B9) 

This gives the result 



dx 



1 







2(27r)3 



(BIO) 



as discussed in Sec. IV Bl 

The difference in support of the light-cone and specta- 
tor integrals is due to the different way the singularities 
of the triangle diagram are handled. An analysis of the 
advantages and disadvantages of these two methods is 
beyond the scope of this paper and will be discussed in a 
future work. Here note only that from the standpoint of 
the spectator method, the integral over kz is limited by 
the light-cone requirement that < a;, which translates 
into the limit < Eg — kz, implying 



m, 



+ fcj - 
2M 



< kz 



This limit breaks the symmetry between the components 
of k, with kz treated differently from the other two com- 
ponents. If one is not careful, the light-cone method 
breaks rotational invariance, and can cause problems at 
low energy where manifest rotational invariance is an im- 
portant constraint. 



APPENDIX C: DETAILS OF THE 
CONSTRUCTION OF THE TWO-PION CUT 
TERM, EQ. (ffSj) 

Unfortunately, chiral calculations of the nucleon form 
factors [46l.l47j only give information about the imaginary 
part of the isovector form factors in the timelike region, 
over a small range of near the onset of the two-pion 



cut (in the region of 4to^ < t ^ q 



-Q2 < 



The real parts of the form factors for both spacelike 
and timelike regions are obtained using dispersion the- 
ory, and are sensitive to the high momentum behavior 
(t = — Q2 j^^2-| q£ ^j^g imaginary part. 

Hence, to use the results of the chiral calculations we 
chose a mathematical function that reproduced the chiral 
calculations of the imaginary part at small t, but gave a 
real part with a strength that could be varied. A con- 
venient function with a simple behavior in the complex 
plane is 



9.-W = 



''.(l-^)(l-^)--.4i 

with the width function 

lit) 



2 3 



(CI) 



(C2) 



where /i = m-^ and — ^tj [j?- — 4 is the relative mo- 
mentum of the two pions (in units of /i) in the tttt rest 
system. Note that the width goes as g^^, as required by 
the P-wave nature of the p, and that (7i_(0) = 0, so 
does not contribute to the charge or magnetic moments. 
The imaginary part of gi_ for t > 4/1^ is 



Imgi_(t) = 



7(t)/t 

' m 



hj t^E{t) 



where 



^ 1 - -2 1-32 



(C3) 



(C4) 



<< 



(Bii) e 



The expansion holds for low (/tttt momentum (q^r 
M), when E{t) ~ 1. This shows that Imgi_(f) is inde- 
pendent of hi to lowest order, as shown in in Figs. [7] and 
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As it turns out, the simple function "/{t)/t does an 
excellent job of fitting the chiral calculations, provided 
we choose qe = 0.17 and aM = 0.68. The mass — 
28y^^ was chosen to equal the square of the p mass, and 
771-2 = 50/i^ is large enough to provide the needed high 
momentum convergence without seriously affecting the 
shape of the imaginary part in the low t region. 



Finally, (jCip may be analytically continued to t = 
— < 0. Giving t a small positive imaginary part and 
using 

i^{t + ie)/Ai2 _ 4 V4 + QVm' (C5) 
gives the result ([75)) . 
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